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I. INTRODUCTION 

Electron and positron fluxes in Cosmic Rays (CRs) have been measured by many 
experiments. For energies above about 10 GeV the positron fraction displays a rising 
behavior. The total flux also displays some features around 100 GeV. There is a common 
agreement on the fact that these data cannot be interpreted solely in terms of known 
astrophysical sources. An unknown source of electrons and positrons CRs has to be 
introduced in order to account for the excesses seen on the top of the background fluxes, 
associated to known astrophysical sources. We refer to [^ E] for recent reviews. 

In [3] we suggested a strategy to investigate the possible charge asymmetry in the 
electron and positron CRs excesses; as a result, we showed that, at that time, even large 
deviations from charge symmetry were experimentally viable. This kind of analysis de- 
serves now to be updated in the light of the recent data collected by AMS-02 |1] and by 
Fermi-LAT [5l |6]. This is precisely the goal of this paper. Future experimental observa- 
tions could even better constrain the amount of charge asymmetry. 

The amount of the charge asymmetry of the CR lepton excesses is crucial to understand 
the physical properties of the unknown source. Among the many candidates suggested, 
there are for instance: astrophysical sources, like supernovae or pulsars (see e.g. [1112113" 
[To]), which are expected to be charge symmetric; dark matter (DM) annihilations and/or 
decays. Contrary to the case of annihilating DM, decaying DM can lead to a charge 
asymmetry provided that both charge conjugation and lepton flavor are violated [H]. 

The paper is organized as follows. In sec. [IT] we review the experimental data to be 



used in the analysis. Section III introduces the notations and a useful relation [3]. After 



having discussed background models in sec. IV, in sec. Mwe investigate the experimental 



status of the charge asymmetry. Here we find a hint of an asymmetric excess, favouring 



the electron component over the positron one. We conclude in sec. VI, where we offer 
the physical applications and interpretations of our results. 



II. DATA ON ELECTRON AND POSITRON CRS FLUXES 

We now briefly summarize the experimental data on electron and positron CRs fluxes, 
(j)e-{E) and (j)e+{E), where E is the energy of the detected e^. 

In 2009, the PAMELA experiment yjj measured the positron fraction, 
(f)e+iE)/{(pe-{E) + (f)e+{E)), betwccu 1 and 100 GeV, finding that it unexpectedly in- 
creases above 10 GeV. This rising behavior is difficult to explain via secondary production 
of positrons in interactions of high energy hadronic CRs. Therefore it has been interpreted 
as a positron anomalous excess in the CR energy spectrum above 10 GeV. This would 
imply the existence of an unknown source of CR positrons - see for instance the nice dis- 
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cussion in [8] . On the other hand PAMELA did not observe any excess in the anti-protons 
spectrum [T3] . 

Aheady in 2008 ATIC ^A\ and PPB-BETS [15] reported an unexpected structure in 
(j)e-{E) + ipe+i^E), in the energy range between 100 GeV and 1 TeV. The picture was soon 
corroborated via the higher-statistics measurements by Fermi-LAT [16] and H.E.S.S. |17j . 
that suggested a possible small additional unknown component in the total flux, on the 
top of the standard astrophysical model predictions. The latter generically assume a 
single-power-law injection spectrum of e^. Fermi-LAT [5] determined that the total e^ 
spectrum in the energy range 7 GeV< E < 1 TeV is indeed compatible with a power-law 
of index —3.08 ± 0.05, but it also displays significant evidence of a spectral hardening 
above 100 GeV. 

In 2011, new experimental informations were added. PAMELA measured the electron 
spectrum between 1 and 625 GeV [18]. Fermi-LAT [6] measured the separate cosmic-ray 
electron and positron spectra, thought with a worse sensitivity than the total spectrum. 

Very recently, AMS-02 [3] measured the CR positron fraction with unprecedented 
precision and up to energies of about 350 GeV. This experimental information, together 
with the precise determination of the total flux by Fermi-LAT [5J, at present, calls for a 
new study of the charge asymmetry in the electron and positron excesses. This is what 
we will perform next following the strategy proposed in [3]. 



III. NOTATION AND SUM RULE 

The observed flux of electrons and positrons can be written as the sum of two contri- 
butions: a background component (j)^±{E), describing all known astrophysical sources; an 
unknown component (f)^±{E) (of whatever origin), which is needed to explain the features 
in the spectra observed by experiments. Explicitly, 

0e+ (E) = 0^+ (E) + 0f+ (E) , 0,- (E) = <I>1 (E) + 0f- (E). (1) 

AMS-02 and Fermi-LAT measure respectively the positron fraction and the total electron 
and positron fluxes as a function of the energy E: 

F^E) = tf[^} , T{E) = ct>AE) + <PAE) . (2) 

(pe+{E) +(pe-{E) 

The left-hand side of the equations above refer to the experimental measures. Given such 
data, our aim is to investigate the unknown contribution leading to the lepton excesses: 

0,^+(E) = F+(i5;)r(E)-0f4E), (3) 

<PUe) =T{E) {I ~ F^{E)) - <pl{E) . 



Clearly, this can be done only by assuming an astrophysical background model, as dis- 
cussed below. 

Here we are interested in particular in a fundamental property of the unknown contri- 
bution: its charge asymmetry [3]. The ratio of the unknown electron and positron fluxes 
is a direct measure of such charge asymmetry: 

'^^ ^ - e(E) F^{E) T{E) - e(^) ■ ^ ' 

Note that this equation can be rewritten as a sum rule [3], 

T{E) l-{l + ru{E))F^{E) ^ 



-(^) 1-ME)^^ 



(5) 



<I>"{E) 



that links together the experimental results, the model of the backgrounds and the de- 
pendence on the energy of the charge asymmetry of the unknown excesses. We use the 
E >25 GeV data bins, since the lower energy bins are affected by the solar modulation. 



IV. ASTROPHYSICAL BACKGROUND MODELS 

Primary electrons can come from galactic CRs while interactions of CRs with the inter- 
stellar medium sources secondary electrons, positrons and antiprotons. The propagation 
of the signal and background fluxes from their production region to the detector is affected 
mainly by diffusion and energy losses. We evaluate the background fluxes at Earth using 
the studies [3 fT9V[2T] . These fluxes can be conveniently described by a power low, with 
a global normalization and a spectral index. We discuss below the impact of the spectral 
index uncertainties. 

We model the background spectrum using 



iE) = N^^B,4E) , (6) 



iB,^,_.,B 



where iV^ are normalization coefficients and Bf.±{E) are provided using specific astro- 
physical models. In this paper we adopt various models, in order to study how much the 
results are affected by the background model choice. In particular, measuring E in GeV 
and the -B's in units of GeV~^cm~^sec~^sr~^, we consider the following models. 

• Moskalenko and Strong (MS) [22], a popular model (used also in [231 El]) for which 

4.5E°-^ 
^'^^^^ " 1 + 650^2-3 + 1500^4-2 ' ^^^ 

O.IGE-^-^ Q.7QE^-^ 

Be-[E) — - , 1 1 T-iflQ , o OTT''? 1.^ "•" 



1 + 11^0-9 + 3.2^2.15 I + 110^1-5 + 600^2-9 + 580^4.2 



• Spectral Index (SI) model, a generic model that we parametrize as: 

B,+ {E) = 1.40X10-M — 
B,-{E) = 5A3xlO-''(^ 

where Eo = 33.35 GeV. The normalization coefficients of the SI model have been 
chosen so that the Be± functions equate those of the MS model a.t E = Eq. Note 
also that the SI model with 7e- = 3.21 and 7e+ = 3.41 actually corresponds to 
the Fermi Collaboration (FC) model zero [251 126]. The SI model with 7e+ = 3.5 
approximates very well the positron background of the MS model. The electron 
background of the MS model has spectral index 3.25 for energies above 100 GeV. 

The background models are illustrated in fig. u\ assuming for definiteness N^± = 0.73. 
The natural range of the electron spectral index of the SI model is 7e- = [3.18,3.26], 
while the natural range of the positron spectral index is 7e+ = [3.4, 3.5] (see for instance 
[71 Uni EOj). The FC model can thus be seen as an SI model with intermediate values 
of the spectral indexes. For comparison, the plot also shows the Fermi-LAT data points 
for the total [5] and separate electron and positron fluxes [B]. To make a comparison 
with other recent studies, we note that refs. [23l |2l] focus on the MS model, allowing 
for variations of its spectral index of about 0.05; this corresponds to our SI model with 
7e- = [3.19,3.29] and 7e+ = [3.45,3.55]. 

The energy dependence of the background fluxes is encoded in the Bf.± {E) functions 
but, as can be seen from eq. (|6]), there is also the problem of fixing the normalization 
coefficients iVj, . As we are going to discuss, in order to characterize a certain background 
model, it is not necessary to make two independent assumptions on N^_ and N^+, but 
just one assumption on their ratio: 



K 

m 



rB = ^. (9) 



If we focus on a certain energy bin E where both the positron fraction F^ and the total 
flux T have been measured, eq. (pi) can indeed be used to derive a range for A^JL : 



jyB(^ (E)r )- ^(^^ l-{l + ru{E))F4E) 

^ ' ^ TB B^.{E) 

For a certain background model (characterized by the Bf.±{E) functions) and using the 
experimental data on T{E) and Fj^{E), the allowed range for N^^ can be calculated by 
making assumptions on the values of ru{E) and r^. For the sake of our analysis, we 
consider it safe to let r^ vary in the range [0.5, 2] (as done also in [231 El])- 
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FIG. 1: Models for background fluxes: MS (dashed), SI for extreme (dot-dashed) and inter- 
mediate (solid) values of the spectral indexes. The latter corresponds to the FC model. The 
fluxes have been normalized by choosing A^^ = 0.73. The Fermi-LAT experimental data on the 
the total flux |5j (red) and separate electrons and positrons fluxes [6j (green) are also shown for 
comparison. 



We consider in particular E = 33.35 GeV (for which B^-{E)/ B^+{E) = 38.82) and 
display the results for N^_ in fig. ^ showing the dependence on r^ for fixed values of ru{E) 
in the left panel, viceversa in the right panel. The thickness of the curves is obtained by 
considering the variation of N^_ associated to the la ranges of T{E) and F^{E). The 
variation due to T{E) turns out to be the dominant one. 



V. CONSTRAINING THE CHARGE ASYMMETRY 



Having derived a range for the normalization of the electron background A^^ as a 
function of ru{E) and tb (so that both T and F+ values are reproduced at a certain 
energy E), we can use it to extrapolate the positron fraction at any energy value. We can 
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FIG. 2: Values of A'^ _ according to eq. (10) assuming E = 33.35 GeV. Left: as a function of 



tb for different values of ru{E). Right: as a function of ru{E) for different values of tb- The 
thickness of the bands is obtained by considering the la ranges of F+(E) from AMS-02 [4j and 
T{E) from Fermi-LAT [5]. 



in fact rewrite eq. ^ as follows: 



F+{E) 



l + ME) 



l~N^.{ru{E),rB] 



Be-{E) 
T{E) 



1- 



ru{E)B,+ {E) 
TB B,^{E) 



:iii 



Suppose now that we specify a background model (namely B^+{E),Be-{E),rB) and that 
we make an assumption on ru{E): it is then possible to check the consistency between 



the extrapolation for different energy values based on eq. (11 ) and the experimental data. 
Clearly, if we use the range of values for N^.{ru{E),rB) discussed previously, we are 
guaranteed that both the total flux and the positron fraction reproduce the experimental 
data in the E energy bin. 



A. Energy independent rij(E) 



In order to test whether current CRs data could support charge asymmetric lepton 
excesses, as done in [3], the first step is to consider the oversimplifying assumption that 
ru is nearly constant in the energy region of interest. 

The extrapolation of the positron fraction, assuming that ru remains constant over 
the entire energy range (from about 30 GeV up to about 700 GeV), is shown in fig. [3] for 
the FC (shaded) and MS (dashed) background models and by taking r^ = 1. We focus 
in particular on the cases rjj = 0, 1/2, 1,2,4 (from top to bottom). The thickness of the 
curves correspond to the la ranges of AMS-02 |1] and Fermi-LAT [5] . One should not be 
worried by the wiggles between in the 30 — 50 GeV, as they are simply due to the features 
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of the data points in that range. The picture shows that the current data are consistent 
with r;7 = 1 for energies up to 100 GeV, but favor a deviation from charge symmetry 
at energies above 100 GeV, where the preferred charge asymmetry value is rj/ ~ 2: the 
latter value means that the electron excess of unknown origin should be about twice the 
positron one. 
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FIG. 3: Positron fraction F^{E) for different values of ru, according to eq. (11). The shaded 
(dashed) curves refer to the FC (MS) background model with rg = 1. The AMS-02 data [1| with 
la error bars (statistical and systematic combined in quadrature) are shown for comparison. 



These interesting results makes it mandatory a deeper study of their dependence on 
the model background: this can be done by considering the impact of varying r^ and the 
spectral indexes 7e±. 

Focusing on the FC background model for definiteness, the top panel of fig. |4] displays 
the dependence on rs- Lowering tb goes in the direction of reducing the positron fraction, 
alleviating the tension between the AMS-02 data ^ and the charge symmetric case above 
100 GeV. However, the global shift turns out not to be strong enough to account for 
ru = 1- 

As a second study of the robustness of the deviation from charge symmetry above 100 
GeV, we consider a generic SI model with r^ = 1 and analyze its dependence on j^- and 
7e+ : this is done respectively in the middle and bottom panels of fig. |4j We can see that 
while 7e- has a significative impact on the slope of the positron fraction, ■y^+ does not 
affect it much. The tension with charge symmetry above 100 GeV is nearly removed for 
values of 7e- smaller than 3.18. 

As a further test of our results, we perform a chi-squared test using the following test 
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FIG. 4: Positron fraction F^{E) for different values of rjj. The AMS-02 data [1] with la error 
bars are shown for comparison. Top: dependence on rs for FC model. Medium: dependence on 
7e- for SI model with 7e+ = 3.41 and r^ = 1. Bottom: dependence on 7e+ for SI model with 
7g- = 3.21 and rs = 1- 



function. For the background we assume the following SI form, 

(/.f±(E)=iVf±E-^^±, (12) 

while for the unknown source we consider 

0^± (E) = N^± E-^^ e"^/^^ . (13) 

The associated charge asymmetry is constant and parametrized by rjj = N^_/N^+. We 
fit the AMS-02 [Ij and Fermi-LAT data [5J to the following seven free parameters for a 
given ru'- 

jfB~, 7e--7e+, jfs^ lu-le-, Eu, N^^ , 7^- , free parameters . (14) 

A reasonable agreement with data, according to the chi-squared distribution, can be 
obtained for any value of ru between 1 and 2. We find that for lower values of ru, 7e- 
turns out to be small, around 3.14 — 3.16 at 90% CL. for ru = 1. For vu = 2 instead the 
7e- range is 3.205 — 3.225. These results are in agreement with our previous comments 
that a value of 7e- closer to the astrophysically expected range 3.18 — 3.26, is obtained 
for Tu ~ 2. 



B. Deriving ru{E) from data 

For various background models, we now study directly the charge asymmetry of the 
unknown excesses, ru{E) = (j)^[E) / (j)^[E) , by considering its expression given in eq. (H), 
which we report here: 

.p. T{E) (1 - F^{E)) - N^.{ru{E),rB) B,-{E) 
F+{E)T{E)- ''^-^Z7 BAE) 

Such study can be done by assuming a background model (namely B^+{E) , B^-{E) ,r b) 
and making an assumption on ru{E). 

In the top panel of fig. [SJwe display the ru{E) range for the MS (dashed), FC (shaded) 
and SI (dot-dashed) models, assuming r^ = 1 and ru{E) = 1. For the SI model we 
choose 7e- = 3.18 and ■je+ = 3.41. The thickness of the curves corresponds to the la 
ranges of AMS-02 |1] and Fermi-LAT [S] . Despite the oscillations below 50 GeV (which 
are simply due to the pattern of the experimental data points), one can see that ru{E) 
displays an increasing behavior with energy. A transition occurs above 100 GeV, where 
ru{E) becomes bigger than unity, spanning the range between 1 and 2, for the MS and 
FC models. The SI model with a low value of the electron spectral index, 7e- = 3.18, is 
instead compatible with charge symmetry. 

9 







dashed: MS, shaded: FC 




i.i 


dot-dashed: SI with 7^+=3.41,y^~=3. 18 /'~\ 




?,n 


: r,=l 






^ — -^ ^^ ^- 








'^ ~' / \ 




, ^ 




; / 




w 


1.5 


/ -^ y 




C' 




- rij{E)=\ / --' /'--^-_ .--"' ^^ 






1.0 


^ '-' - '' 






f\ '~^. ' ,c-^-^y ^-^ ------- 






- II \ /Cv^ --^' ~' 1 "-—'" 




0.5 
n n 





50 



100 
E [GeV] 



200 



300 



4- 






- 


\ 










^\ 




- 


/ 


//"" ' " -";;■ _ - - - 




/ ''\.4'^{EpfdS"'"-'' 


FC withris=l 



50 



100 
E [GeV] 



200 



300 



FIG. 5: Ratio ru{E) according to eq. (15). Top: for the MS (dashed), FC (shaded) and SI 
(dot-dashed) models with r^ = 1 and ru{E) = 1. Bottom: for the FC model with r^ = 1 and 
various values of ru{E). The thickness of the curves corresponds to the la ranges of AMS-02 
^ andFermi-LAT [5]. 



Since we are particularly interested in deviations from charge symmetry, it is important 
to study the dependence of ru[E) on the value chosen for ru{E). This is done in the 
bottom panel of fig. [5} focusing on the FC model with r^ = 1. One can see that for 
whole interval ru{E) = [0.5,4], ru{E) is bigger the unity above 100 GeV. In addition, 
rjj{E) ~ 2 is the sole case that allows rjj{E) to be nearly energy independent. For 
ru{E) > 2 {rjj{E) < 2), ru{E) is a decreasing (increasing) function of the energy. We 
also considered the dependence of the ru{E) on the value chosen for r^, finding it to be 
negligible. 
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VI. APPLICATIONS AND CONCLUSIONS 

The physical sources of the excesses observed in the electron and positron CR fluxes 
are naturally divided in two classes: charge symmetric and charge asymmetric. The first 
class is characterized by rjj = 1. All other sources, for which rjj{E) depends on the energy 
or is constant but different from unity, belong to the second class. 

This classification applies indiscriminately to any source, be it astrophysical or of DM 
nature. It is not the goal of this paper to enter in the merit of any specific model. Here we 
merely classify the main models suggested in the literature with respect to their potential 
to yield a charge asymmetry. Then we suggest how to use our results about charge 
asymmetry, expressed in fig. |5| for a straightforward test of any production mechanism. 

Astrophysical models 

The simplest models of pulsars - see e.g. [9] and references therein - are charge symmetric, 
given that the basic assumption is that pulsars inject the same number of electrons and 



positrons in the interstellar medium. Our parameterization eq. (13) with N^+ = N^_ 
(hence ru = 1) describes well the fluxes from a pulsar expected at Earth. 

For supernovae [21 [71 El [10] the situation is more delicate. There is a primary source 
of electrons which is responsible for Bf.-{E) and, on top of that, an equal number of 
positron and electron secondaries produced at the source and which might be responsible 
for the excesses. Here too the production mechanism for the excesses should correspond 

to rjj = 1. 

DM models 

Symmetric DM can lead to CRs by either annihilation, decay or both. However, as proven 
in [TT], all models of symmetric DM imply rjj = 1 for any energy. A recent re-analysis 
has been performed in ref. [23]. 

In order to achieve an ru ^ 1, DM must be asymmetric (therefore decaying) [271 [28] 
and furthermore violate lepton flavor symmetry [TT1[25]. These models lead to an energy 
dependent ru. For instance, for an asymmetric DM candidate decaying into yU^r+ we 
obtain a naturally increasing behavior for ru{E), from about 1.5 at -E = 30 GeV up to a 
value of 3 at E = 300 GeV, see fig. 3 of [IT]. 

For the DM interpretation of the excesses, however, attention must be paid to gamma- 
ray constraints, see for instance 



We have shown that by combining the recent data from AMS-02 with those from 
Fermi-LAT, a charge asymmetry in the unknown excesses of electron and positron CRs is 
emerging, favoring the electron component. The result relies on having chosen a conser- 
vative estimate for the astrophysical background fluxes. Charge symmetry can be rescued 
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when adopting an electron background spectral index slightly smaller than the commonly 
assumed values. 

Given that both astrophysical and DM models can be classified according to their 
predictions for the asymmetry, the impact of the future AMS-02 [3] results for the charge 
asymmetry will play a crucial role in discriminating the proposed production mechanisms. 
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